Landau problem for bilayer graphene: 
Exact results 
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We consider graphene bilayer in a constant magnetic field of arbitrary orientation (i.e. tilted with 
respect to the graphene plane). In the low energy approximation to tight binding model with Peierls 
substitution, we find the exact spectrum of Landau levels. This spectrum is two-fold degenerate in 
the limit of large in-plane field, which gives rise to a new SU(2) symmetry in this limit. 
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I. INTRODUCTION 

There are not many examples of exactly solvable prob- 
lems in physics. It is even more intriguing, when such a 
problem emerges in a hot topic with a rich perspective 
of potential applications. Graphene, one or a few layers 
of two-dimensional hexagonal carbon lattice of graphite 
[l[ , attracted enormous interest in recent years due to its 
unusual properties, like (quasi)relativistic quasiparticle 
spectrum as well as prospective applications to construc- 
tion of electronic devices (see |2j for a recent review). 
Among these, the bilayer graphene is believed to be one 
of the most rich in features at the same time accessi- 
ble technologically. In particular, it exhibits a gap in 
the electronic spectrum, which can be tuned by exter- 
nal electromagnetic fields. This explains why the bilayer 
graphene in external fields, in particular, in a magnetic 
field is a subject of intensive recent study @-|(|. 

In the present letter we consider the general Landau 
problem for the bilayer graphene system in a constant 
magnetic field having an arbitrary direction with respect 
to the graphene plane. We show that the energy spec- 
trum can be found exactly in terms of Spheroidal func- 
tions, which are solutions to a particular case of Con- 
fluent Heun's equation 0], and corresponding eigenval- 
ues. Further analysis of the eigenvalues reveals that in 
the limit of strong in-plane magnetic field the spectrum 
becomes two-fold degenerate, which leads to additional 
SU(2) symmetry. 

The plan of the remainder of this note is as follows. 
In the next section we introduce, starting from the tight- 
binding Hamiltonian, the low energy Lagrangian describ- 
ing the electronic wave function in bilayer graphene in 
external electromagnetic field. Then, we consider con- 
stant magnetic field. As a warm-up exercise, we derive 
the Landau level (LL) spectrum in perpendicular mag- 
netic field. After this we consider the generic case of 
tilted magnetic field, and find the eigenvalue spectrum. 



Finally, we analyze the spectrum and reveal an asymp- 
totic two-fold degeneracy as well as compare it to the 
purely in-plane magnetic field result. The technical de- 
tails behind the results reported here will appear in an 
accompanying paper Q . 



II. BILAYER GRAPHENE IN CONSTANT 
MAGNETIC FIELD 

According to the tight-binding model the bilayer 
graphene with Bernal stacking can be described by the 
following Hamiltonian, 
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where a, b and their conjugate ST, w are, respectively an- 
nihilation and creation operators for A and B sites of 
the upper layer, while lower layer operators carry no 
tildes. The indices n run through A/B lattice sites, S a , 
a = 1,2,3 are the vectors connecting this site to the near- 
est neighbors on the upper layer, — d a does the same for 
the lower layer and a — ± is the electron's spin. Param- 
eter 70 is the nearest neighbor hopping amplitude, and 
71 is the leading interlayer tunneling amplitude (see fig. 
[!}. The experimental values of the couplings found e.g. 
in are, 
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2.9 eV, 
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0.3 eV. 



(2) 



Electromagnetic interaction is turned on through the 
Peierls substitution, i.e. any two operators located at 
different sites, e.g. n and n + 6 are connected by a phase 
factor representing the parallel transport in electromag- 
netic field, 
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FIG. 1: Schematic description of the graphene bilayer lattice 
with Bernal stacking, a cc w 1.42A and d « 3.35A. The dotted 
lines represent the interlayer tunneling. We use tildes for the 
parameters referred to the upper layer and no tildes for the 
lower layer ones. 



If the in-plane component -By vanishes, i.e. the mag- 
netic field is perpendicular to the graphene plane, B± = 
B, the Hamiltonian §5§ reduces to, 
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where m = 71 /2V 2 , ps 0.030 x m c is effective mass of the 
fcrmion. 

The operator H can be readily diagonalized by stan- 
dard techniques. Indeed, introduce the oscillator rais- 
ing/lowering operators a and ot* with standard commu- 
tation relations according to 



At low energy the quasiparticle is described by an eight 
component wave function \&(x, t): There are two compo- 
nents for each Dirac valley ±K as well as two values of 
the spin index. Although a subject to quadratic disper- 
sion relations, the quasiparticle is a chiral particle. 

In the case of constant magnetic field the effective de- 
scription can be casted into the following Lagrangian, 
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where the chiral 'Hamiltonians' Wr±jg\ correspond to dif- 
ferent Dirac points ±K, where K = (47r/3\/3acc, 0) f]~3j j . 
The operators T-L(+tC\ are most conveniently expressed 
using the complex coordinate for the graphene plane, 
w = -^=(x + iy). Thus %(k) is given by 
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while the operator "H(_k) is obtained from %(k) by inter- 
change of D and D. In eq. ((5|) we neglected the Zeeman 
term, the only effect of which is a shift of energy eigenval- 
ues by ±(he/2m e c)B. The covariant derivative D and its 
conjugate D are defined, respectively, as D = d w + j^A w , 
D = dw + ^A»i A w and being the complex compo- 
nents of the vector potential. Peierls' phase 
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corresponds to the parallel transport between the layers. 

In the case of a constant magnetic field we can fix the 
gauge such that electromagnetic potential A, is given by 

A w = A% = -^B ± w, A z = -J=S|| [r)w - fjw). (7) 

Here B± and B^ij are, respectively, perpendicular and in- 
plane components of the magnetic field, the phase factor 
•q gives the direction of the in-plane projection of mag- 
netic field and can be absorbed into redefinition of w. 
We restrict our analysis to a single Dirac point K and, 
respectively, drop the subscript. 
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Then the problem is equivalent to the eigenvalue problem 
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where s corresponds to the energy eigenvalue E as E = 
— (HeB/mc)e. 

The solution to the eigenvalue problem (|10[) is given 

by 



E„ -- ±h [ — ) y/n(n - 1), 



(11) 



parameterized by n — 0,1,2,... and the sign ±. Let 
us note that due to chirality the Hamiltonian is not sign 
definite. 

Now let us turn on the in-plane magnetic field and con- 
sider the modifications in the eigenvalue problem. The 
Hamiltonian now can be written as 
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where v = (edB^/Hcy/2). One can easily find the zero 
modes of this operator, so let us concentrate on the 
nonzero eigenvalues. In this case the eigenvalue prob- 
lem for Hamiltonian (fl~2"|) can be reduced to the following 
one component problem, 



a f (a f +P)(a + 0)oeipi = Xtpi, 



(13) 



with /3 = Vyl an( l ^ = £2 ■ The remaining component 
ip2 is expressed through, 

^a = A- 1 e"( ffl - ,0 )(a + j9)Q^i. (14) 

The one-component problem (|13l) is not as simple as in 
the purely perpendicular case, so let us consider it more 
closely. Using the antiholomorphic representation: 
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FIG. 2: A plot of Energy eigenvalues (eV) as a function of in- 
plane magnetic field (T) for the fixed value of perpendicular 
field B± = 10~ 4 T. Red lines represent even energy levels, 
while blue lines the odd ones. One can clearly see convergence 
of even and odd neighbor eigenvalues. Zeroth and first states 
are the zero modes of the Hamiltonian, therefore, they are 
degenerate at any value of field. 

the eigenvalue problem is reduced to the second order 
differential equation 

x(x + p)(u" + pu') - Xu = 0, (16) 

which we can solve on the real line and then analytically 
continue. The finite norm condition for eigenstate ip im- 
plies that function u should be a regular function in the 
interval x € (—00,00) with a moderate growth (at most 
exponential) as x — > ±00. 

The equation (1 16[) is a particular case of confluent 



Heun's or spheroidal equation. By a change of indepen- 
dent variable x n- {P/2)(x— 1) and then of the dependent 
one u(x) = y/x 2 — 1 exp(— ^-x)v(x) it reduces to one of 
the standard forms of the spheroidal equation. The solu- 
tions satisfying physical conditions are given by the an- 
gular oblate spheroidal functions, 

where c n (/3) is the normalization constant. These func- 
tions represent Landau eigenstates. The correspon ding 
Landau levels are given by the spheroidal eigenvalues [14| 

where the superscript (a) stands for 'angular'. Recall, 
that P = y/ed 2 /2hc (By/VST). 

Spheroidal equations appear in many situations, like 
molecular hydrogen ion eigenvalue problem, radio an- 
tenna description as well as geophysical applications, 
therefore spheroidal functions and respective eigenvalues 
are relatively well studied [?J . They are also implemented 
as standard libraries in Mathematica [101 ] . 

A typical behavior of energy eigenvalues is depicted in 
Fig. [5] From it one can see, that for f3 = the eigen- 
value spectrum is reproducing that of strictly perpendic- 
ular magnetic field problem. These eigenvalues are quite 
stable against small in-plane field perturbations. Indeed, 
small P expansion for the spheroidal eigenvalues gives, 
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128(2n + 3)(2n + l) 3 (2n - 3) 3 (2n - 5) 



/3« + 0(/3 iu ), (19) 



so, the leading correction is of the order ~ S| / B\ , which 
means that for a reasonably small tilt the LL spectrum 
is very close to the purely perpendicular magnetic field 
values. 

In the opposite case, when the magnetic field is almost 
parallel to the graphene plane, we can use the large /3 
asymptotic expansion of the eigenvalues [ll| , 
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where the square bracket denotes the integer part. The 
presence of this function implies that there is an asymp- 
totic degeneracy between the each even level and the fol- 
lowing odd one. The difference between these levels is 



exponentially vanishing: 
2 
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where k — 0, 1, 2, . . . This degeneracy can be interpreted 
as emergence of a new SU(2) symmetry in the nearly 
in-plane limit of the orientation of magnetic field. This 
asymptotic symmetry is mixing the neighbor states with 
numbers 2k and 2k + 1. 

Although, mathematically even and odd eigenval- 
ues are quickly converging to each other, in practice, 
the effect can be caught only at very large values of 
the magnetic field and/or very precise in-plane orien- 
tation. Indeed, for magnetic field measured in tesla, 
P ~ 10~ 3 i3|| j\JB^. Assuming that degeneracy for the 
lowest levels occurs at the values P ~ 10, we have 
89 ~ 10 -7 £>, where 59 is the angle deviation from the 
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parallel direction. Currently available magnetic fields are 
of the order ~ 10 2 T, which means that to capture a mea- 
surable effect the magnetic field should be parallel to the 
graphene plane with the accuracy ~ 10 _5 rad. 

More generally, the asymptotic expansion formula (|20)) 
predicts that in the limit in which in-plane component is 
much larger than the perpendicular one the energy eigen- 
values scale like E n ~ Bu y/B± or in terms of magnitude 

and angle as ~ B 3 / 2 y/\6\. Thus, the energy eigenvalues 
in such a magnetic field exhibit a hybrid behavior reveal- 
ing both relativistic and non relativistic patterns. Thus, 
if we fix the (large) value of the in-plane Bt\ the spectrum 
scales according to the relativistic rule: E ~ \/B±. 

It is instructive to compare the above limit with the 
exact solution in parallel magnetic field, which we can 
obtain after setting B^ = 0. 

In the parallel magnetic field the Dirac point Hamilto- 
nian ([5]) becomes, 

V __K* ( 0e-"(*-™>0\ (99] 

This operator can be easily diagonalized using the plane 
wave basis uj, k ~ e -(kw-kw)^ wn i cn gives a continuous 
spectrum of eigenvalues parameterized by the wave num- 
ber k: 

El, k = ^\k\ 2 W-k\ 2 (23) 

where v is defined below the equation (|I2[) . In addi- 
tion to the invariance with respect to inversions about 
x-axis, the spectrum is left invariant by the transforma- 
tion k H> v — k. In particular, there are two zeroes of the 



dispersion relation (|23|) : at k = and k — v. The group 
of symmetry transformations mixing the waves u k k and 
u v-k.v-k is SU(2). 

In the low energy regime (in this case, when E <C 
^-^-v 1 ) the energy can be expanded around the zeroes 
of dispersion relations: k = and k — v. This yields two 
massless Dirac fermions separated by momentum shift 
v as discussed in [l2j]. The quantum numbers counting 
these 'vacua' are related to the SU(2) group which de- 
scribes the symmetry of the spectrum. In this limit the 
whole electronic system is described by a sixteen compo- 
nent wave function. 

It is clear, that this SU(2) symmetry can be identi- 
fied with the asymptotic SU(2) in the large /? limit of 
the tilted magnetic field system. The easiest way to con- 
vince ourselves in this is to consider the zero eigenvalue 
levels of the tilted problem. Unlike the higher levels, 
they can be expressed in terms of elementary functions. 
There is a simple relation between the zeroth eigenvalue 
eigenfunctions ipo(w,w) and ipi(w,w) — e^ Qt ipo(w, w) — 

y he 

e " B ± iftoyw, w). This relation implies that, in the par- 
allel field limit B± — > 0, the Fourier modes are related 
up to a constant phase factor by the transformation: 
k i y v k. 
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